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Abstract. While the propagation of tsunamis is weh understood and weh simulated by 
numerical models, there are still a number of unanswered questions related to the gener- 
ation of tsunamis or the subsequent inundation. We review some of the basic generation 
mechanisms as well as their simulation. In particular, we present a simple and computa- 
tionally inexpensive model that describes the seabed displacement during an underwater 
earthquake. This model is based on the finite fault solution for the slip distribution under 
some assumptions on the kinematics of the rupturing process. We also consider an unusual 
source for tsunami generation: the sinking of a cruise ship. Then we review some aspects 
of tsunami run-up. In particular, we explain why the first wave of a tsunami is sometimes 
less devastating than the subsequent waves. A resonance effect can boost the waves that 
come later. We also look at a particular feature of the 11 March 2011 tsunami in Japan - 
the formation of macro-scale vortices - and show that these macro-scale vortices can be 
captured by the nonlinear shallow water equations. 



the water surface is perturbed. There are three main types of disturbances: underwater 
earthquakes concentrated in zones where there is shpping or subduction of tectonic plates, 
submarine landshdes which are often but not always triggered by earthquakes, and sudden 



Contents 



1. Introduction 

2. Tsunami generation 

2.1. Finite fault 

2.2. Costa Concordia 

3. Tsunami inundation 

3.1. Run-up amplification 

3.2. Whirlpool like effect in tsunamis 

4. Conclusion 
Acknowledgements 
References 



1 

2 
4 
5 

14 
14 
14 
16 
18 
18 



1. Introduction 
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earth surface movements adjacent to the ocean (volcanoes, rock falls, sub-aerial landslides, 
ship sinking) . In the generation of tsunamis by earthquakes the key point is to predict the 
displacement of the sea bottom, and then to understand the energy transfer to the water 
column. The generation of tsunamis by submarine landslides is even more challenging for 
various reasons: lack of data, coupling between fluid and solid motions, longer duration, 
more physical parameters. The modeling of sudden earth surface movements adjacent to 
the ocean has not been much studied, the reason being that such movements are quite 
rare. 

The propagation of tsunamis is now well understood and operational codes can easily 
propagate tsunamis across a whole ocean (see for example the various simulations of the 
megatsunamis of 2004 in the Indian Ocean and of 2011 in the Paciflc Ocean). 

Wave run-up (maximum vertical extent of wave uprush on a beach above still water 
level) and wave inundation (maximum horizontal extent) have been studied during the 
last sixty years, but continue to be a challenging problem. A high level of mesh reflnement 
as well as high resolution bathymetric and topographic data are required to describe local 
wave run-up. 

Our group has developed a novel tool for tsunami wave modelling. This tool has the 
potential of being used for operational purposes: indeed, the numerical code VOLNA is 
able to handle the complete life-cycle of a tsunami (generation, propagation and run-up 
along the coast). The algorithm works on unstructured triangular meshes and thus can 
be run in arbitrary complex domains. A detailed description of the flnite volume scheme 
implemented in the code as well as the numerical treatment of the wet/dry transition can 
be found in [7]. 

In Section 2, we review some of the basic generation mechanisms as well as their sim- 
ulation. In particular, we present a simple and computationally inexpensive model that 
describes the seabed displacement during an underwater earthquake. This model is based 
on the flnite fault solution for the slip distribution under some assumptions on the kinemat- 
ics of the rupturing process. We also consider an unusual source for tsunami generation: 
the sinking of a cruise ship. In Section 3, we review some aspects of tsunami run-up. In 
particular, we explain why the flrst wave of a tsunami is sometimes less devastating than 
the subsequent waves. A resonance effect can boost the waves that come later. We also 
look at a particular feature of the 11 March 2011 tsunami in Japan - the formation of 
macro-scale vortices - and show that these macro-scale vortices can be captured by the 
nonlinear shallow water equations. 

2. Tsunami generation 

The modelling of tsunami generation was initiated in the early 1960's by the prominent 
work of Kajiura [14], who proposed the static approach which is still widely used by the 
tsunami wave modelling community: the static sea bed displacement is translated towards 
the free surface as an initial condition. The most classical solution for the co-seismic sea 
bed displacements is the celebrated Okada solution [22]. The kinematics of earthquakes is 
relatively well understood. The maximum bottom deformation is achieved during a flnite 



ON THE MODELLING OF TSUNAMI WAVES 



3 



time known as the rise time. For example the rise time was 8 seconds for the July 17, 2006 
Java event simulated by Dutykh et al. [6]. In the next subsection, we will show how 
kinematics can be taken into account in tsunami generation. 

However, modelling tsunamis generated by landslides is far more complicated. Firstly, 
the kinematics of underwater landslides is not well understood. Also the time scale over 
which they occur is longer than for an earthquake, so simply transferring the sea bed 
deformation directly to the free surface does not accurately model the landslide. For one- 
dimensional landslides, Liu et al. [15] showed that the free-surface elevation is significantly 
diflFerent from the bottom motion by using an anlytical technique developed by TuCK & 
Hv^ANG [35]. Sammarco & Renzi [28] extended the results to two-dimensional landshdes 
(see also [26, 27]). Sarri et al [29] used the Sammarco & Renzi model to build a 
statistical emulator. 

Mass movements on dry land can be put into a number of different categories. Rota- 
tional slides, translational slides, block slides, falls, topples, debris fiows, debris avalanches, 
earthfiows, creep and lateral spreads each have their own characteristic kinematics of mo- 
tion. However, when considering landslides that are in contact with water, our knowledge 
of the kinematics in this environment is lacking. This is due to lack of data, in particular 
lack of bathymetry data prior to large tsunamigenic events. 

A recent review of submarine mass movements [16] shows that they can consist of soil and 
rock, and similar to dry land movements they can take the form of slides, spreads, fiows, 
topples or falls, but in addition they can develop into turbidity currents. They can move 
up to to 50 km/h, reach distances over 1000 km and volumes can be enormous, the largest 
known being the Storegga slide at approximately 2500 km^. Underwater mass movements 
pose a threat to coastal communities and infrastructures both onshore and offshore. The 
main triggers are seismic shaking, overloading gas hydrate dissolution and excess pore 
pressure, wave loading, erosion and human activities such as coastal construction. The 
1998 Papua New Guinea tsunami generated a renewed interest in tsunamigenic landslides. 

Although generally landslide tsunamis occur over a much smaller scale than earthquakes 
(O(lkm) vs 0( 100km)), they can cause very large run-up values. Gutenberg [11] was 
one of the first to suggest that tsunamis can be caused by submarine landslides. He goes 
so far as to say that they can be considered as one of the chief causes of tsunamis. He 
reported that co-seismic landslides possibly triggered large waves in Ceram in 1899, Assam 
in 1897 and caused cable breaks in Greece in the 19th century. More recently, in 1979, a 
part of the Nice harbour extension slumped into the Mediterranean and was followed by a 
small tsunami [1]. 

The Storegga slide is one of the largest known submarine landslides (2500 — 3500 km^) 
and occurred off the west coast of Norway generating a huge tsunami 8200 years ago [2] . In 
1929 an earthquake at the edge of Grand Banks, Canada, triggered a large submarine slope 
failure (200 km^) generating a tsunami with run-up heights up to 13 m that propagated 
as far as Portugal and the Azores Islands [9]. A tsunami generated by a 0.03 km^ of rock 
falling from 914 m into Lituya Bay, Alaska, in 1958 is dubbed the world's largest tsunami 
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at 524 m. Plafker [24] points out numerous major landslides that were triggered during 
the 1964 Alaska Earthquake. 

There is debate about the triggering of the 1946 Aleutian tsunami that destroyed Scotch 
Cap's lighthouse. Freyer et al. [10] argue the possible involvement of a submarine slump 
that caused nearsource damage. Similarly, landslides have been suggested as a possible 
cause of amplification of the 1992 Flores Island tsunami [12]. Ma et a/. [17] demonstrate 
large scale slumping combined with faulting as a good model for the 1975 Kalapana, Hawaii 
tsunami. 

Seafioor mapping has begun to open up the research of submarine landslides. Evidence 
of several underwater slides in the St. Lawrence estuary, Canada, have been identified 
through seafioor mapping [25]. 

The scale of the 1998 Papua New Guinea tsunami was a "wake-up call" for tsunami 
scientists. Although there was much controversy initially over the main trigger of the 
tsunami, the magnitude 7.1 earthquake was relatively small compared to the 10 — 15 
m tsunami that devastated the coast. Synolakis et al. [31] presented high resolution 
bathymetric data combined with hydrodynamic modelling indicating a large underwater 
slump in the area. Subsequent research of this event has contributed to new and improved 
models of tsunamis due to submarine mass failure and it is now accepted by most scientists 
that the source of this event was a submarine slump [32]. 

There is also strong evidence that Island volcanoes such as Stromboli, Italy [33], Ritter 
Island, Papau New Guinea [37] and the Hawaiian Islands [20] have experienced lateral fiank 
collapses in the distant past which potentially could have caused large scale tsunamis. On 
the 30*^ December 2002 part of the western fiank of the Stromboli volcano slid into the 
sea initiating two tsunamis seven minutes apart [16]. This gave TiNTi et al. [34] and 
Chiocci et al. [3] the unique opportunity to observe and model both the submarine 
and subaerial morphological changes thanks to a previous multibeam survey of the area. 
Recently, scientists have become concerned that a section of the Cumbre Vieja volcano. La 
Palma in the Canary Islands, may experience failure and generate a megatsunami [4, 36]. 

Landslides have the potential to cause more destructive tsunamis than originally thought. 
They may be seismic triggered landslides amplifying waves originally generated from earth- 
quakes, large scale collapsing of volcanoes into the sea, or rock falls from mountainous 
regions into restricted bays. In Section 2.2, we will consider a special case of tsunami, 
generated by the sinking of a cruise ship. 

2.1. Finite fault. Here we review some recent advances in seismology and show how to 
reconstruct better co-seismic displacements of a tsunamigenic earthquake. More precisely, 
we use the so-called finite fault solution developed by Ji and his collaborators [13], based on 
static and seismic data inversion. This solution provides multiple fault segments of variable 
local slip, rake angle and several other parameters. By applying Okada's solution to each 
subfault, the sea bed displacement is reconstructed with higher resolution. Since Okada's 
solution consists of relatively simple closed- form analytical expressions, all computations 
can be done efiiciently enough so that they can be used in a real-time Tsunami Warning 
System. Seabed displacements are then coupled with a water wave model. Further details 
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can be found in [6]. This approach was recently extended to include horizontal displace- 
ments as well [5]. Figure 1 shows snapshots of the free surface elevation for the July 17, 
2006 tsunami generation. The x— axis is the longitude while the ^— axis is the latitude. 
The times are in seconds. The water elevation is in meters. 

2.2. Costa Concordia. This is the first application of the VOLNA code presented in this 
paper. A second one will be presented below. The Costa Concordia has been balancing 
nearly horizontally on its side since she partially sank on the night of 13 January 2012 on 
rocks beside the island of Ciglio, Italy (see Figures 2 and 3). The local bathymetry in the 
area contains a steep drop from approximately 20 m depth to over 100 m (see Figures 4 
and 5). Before the ship is stabilised, it is possible that she may slide into deeper waters 
and create waves that would travel along the island and possibly head towards the Italian 
mainland. If this were to happen, could a tsunami be generated? 

A similar incident occurred off the Creek island of Santorini in 2007 with the MS Sea 
Diamond when she ran aground on a volcanic reef. The ship was towed off the rocks and 
stabilised, but the ship sank some 15 hours after it initially struck rocks. The bathymetry 
in the area is a caldera (cauldron-like volcanic feature) and so the shore is almost vertical. 
The stern now sits in about 180 m of water and the bow in about 60m (information obtained 
through Wikipedia). It was feared the wreck would slide deeper into the caldera below. 

We consider here estimates taken from [19], Equation (2.1) and [23], Equation (2.2) for 
wave heights generated by displacements in the bathymetry: 

„,= (sae±'J^Y" (2.1) 

and 

Hlv'^ rexp {—k{x + ct^) — exp {—k{^ + cr)^) 



772 (X,t) = 



2c 



C + V 

exp {—k{x — ctY) — exp {—k{^ — crY)^ 



(2.2) 



c — V 

where ^ = x — vT and r = t — T. In these expressions, r] denotes the free-surface elevation 
(maximum only for rji and space and time distribution for 772)- <^ is the transfer coefficient 
(only a portion of the landslide momentum is imparted on the water column) , Ps is the solid 
density, is the water density, w and i are the two horizontal dimensions of the landslide, 
Hl 18 the height of the landslide, v is the landslide velocity, c is the long wave velocity 
iVdh) with h the water depth, T is the duration of the landslide and A; is a coefficient. 

The dimensions of the ship are 291 x 52 x 38 m^. Assuming that she is lying completely 
on her side, a rough geometry is shown in Figure 6. The following assumptions are used: 
the width parallel to the shore is w = 291 m, the height is = 38 m, the length is ^ = 52 
m, the distance from the top of the ship to the still water line is ^ 38/2 = 19 m, the depth 

is 38/2 = 19 m, and the slope of the shore is ^ 80/400 = 0.2. We also assume that 
Ps = Pw = 1000 kg/m^, a = 0.01 and k = 18/ P. Then the only missing parameter is 
the velocity which of course is difficult to estimate. Some sources say that the Titanic 
sank at approximately 15 m/s once she was completely submerged (information obtained 
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Figure 1. Generation of the July 17, 2006 tsunami in Java. Snapshots of 
the free surface elevation computed with a water wave model. The waves are 
generated by dynamic co-seismic bottom displacements reconstructed using 
the corresponding finite fault solution. 
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Figure 2. Location of the Costa Concordia at time of submission from 
http : //www.marinetraf f ic . com (17/05/12) 




Figure 3. The Costa Concordia uk.reuters.com (retrieved 17/05/12) 

through National Geographic). However, her stern hfted high into the air as the ship 
pivoted down into the water and she may have broken in two before she sank completely. 
Although the Costa Concordia is of the same scale as the Titanic (she is actually bigger) 
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Figure 4. Bathymetric chart from the Itahan Hydrographic Institute. The 
Costa Concordia is near Gabbianara (north of Gigho Porto on the map) 

Costa CQncordia^s precariDus position 




Figure 5. Costa Concordia position from BBC news 
http://www.bbc . co .uk/news/world-europe- 16563562 (retrieved 
17/05/12) 
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Figure 6. Rough dimensions of the Costa Concordia 



the Titanic was in open water and sank bow first, where as the Costa Concordia will have 
friction due to the rocks below her and is positioned lying on her side which may cause her 
to sink differently. 

Applying the above parameters to Equation (2.1) from [19] gives 

(\ 1/2 
(8)(0.01)(l)42i22L 1 = (11.6^)^/^ m. 
^ V(9.8)(19) J 

Sinking speeds of = [0.001, 0.1, 1] m/s will give wave heights of r] = [0.108, 1.08, 3.41] m. 
Therefore, the speed that the ship goes down at is an important parameter to get right. 

Okal & Synolakis [23] model the free surface propagation due to an underwater gauss- 
ian slump source by Equation (2.2). With the above parameters for the Costa Concordia 
together with a total travel time T = 100 s and a velocity v = 1 m/s^ the maximum occurs 
at X = 0, t = 0. The free-surface elevation at this point is 772(0, 0) = 0.2052 m. For = 0.1 
m/s, 772(0, 0) = 0.002 m and v = 10 m/s gives 772(0, 0) = 44.08 m. Again the velocity is an 
important parameter. 

Simulations of the sinking of the Costa Concordia were carried out using the VOLNA 
code. A gaussian shaped slide with the same volume as above (291 x 52 x 38 m^) was 
used as a model of the ship and three simulations were carried out with sinking speeds 
V = (0.5, 1, 10) m/s. The geometry is shown in Figure 7. 

An unstructured triangular mesh was implemented over an area of 8.1 x 16 km^ with a 
distance of 100 m behind the shoreline x = 0. The control volumes were of the order of 
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Figure 7. Geometry used for the VOLNA simulation of the Costa Concor- 
dia on sloping sea floor, from above. 

100 m, but reflned to the order of 5 m in a 600 x 300 m^ region near the ship (see Figure 
8). 

Each simulation was run for 300 seconds. The free surface plots for each simulation are 
shown at t = 75 s, t = 150 s and t = 300 s in Figures 9, 10 and 11 respectively. 

Wave gauges were recorded at three points for each simulation: at the shore where the 
ship sinks (x, y) = (0, 0), on the shore line far from the slide (x, y) = (0, 400) and far offshore 
from the ship at (x, y) = (1000, 0). Figure 12 shows the free surface at (x, y) = (0, 0) from 
the time the ship leaves the shore and is completely submerged. The free surface values 
before this time are not applicable because the water depth is zero. For v = 10 m/s the 
maximum run-up value is 5.2 m, for v = 1 m/s the maximum run-up is 0.52 m and for 

V = 0.5 m/s there is no run-up for t < 300 s. 

Figure 13 shows the wave gauge on the shore line far from the slide at (x, y) = (0, 400). 
Clearly the free surface amplitude for = 10 m/s is much larger than for v = 1 m/s and 

V = 0.5 m/s. The maximum run-up for v = 10 m/s is 2 m, and negligible for the other 
two velocities. 

Figure 14 shows the wave gauge far offshore from the ship at (x,y) = (1000,0). The 
maximum run-up values for ^ = 10, 1 and 0.5 m/s are 3.32,0.12 and 0.05 m respectively. 
Once again these results demonstrate the effect of the sinking velocity. If the ship were to 
sink slowly then a wave of only a few centimeters would be generated, however if it were 
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Figure 8. Mesh used for the VOLNA simulation of the Costa Concordia. 




(a) V — 0.5 m/s 



{h) V = 1 m/s 



(c) V = 10 m/s 



Figure 9. VOLNA simulations of the Costa Concordia sinking: free surface 
at t = 75 s. 



to slide quickly down the steep slope it lies on, waves of a few meters could threaten not 
only the immediate coast, but the Italian mainland and the coastline on the island north 
and south of the ship. In particular, a wave entering the nearby port south of the Costa 
Concordia could be reflected and amplifled. 

Since the sinking velocity is an important factor in determining the size of the waves, it 
is important to know the average density of the ship. Reports on the weight of the Costa 
Concordia vary. She has a dead weight of 10, 000 t but some reports suggest that she has 
an actual weight of 45, 000 t, not including luggage or water inside. If the ship is half flUed 
with water then the mass of the water is 2.88 x 10^ kg and adding this to 45, 000 t would 
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(a) V — 0.5 m/s {h) v = 1 m/s (c) v = 10 m/s 



Figure 10. VOLNA simulations of the Costa Concordia sinking: free sur- 
face at t = 150 s. 




(a) V = 0.5 m/s (b) v = 1 m/s (c) v = 10 m/s 



Figure 11. VOLNA simulations of the Costa Concordia sinking: free sur- 
face at t = 300 s. 
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Figure 12. Wave gauge at {x,y) = (0,0) for VOLNA simulations of the 
Costa Concordia sinking from time when the ship has left the shore. 
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Figure 13. Wave gauge at {x,y) = (0,400) for VOLNA simulations of the 
Costa Concordia sinking: (a) v = 0.5 m/s, v = 1 m/s and = 10 m/s; (b) 
close up of = 0.5 m/s and v = 1 m/s 
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Figure 14. Wave gauge at (x, y) = (1000, 0) for VOLNA simulations of the 
Costa Concordia sinking: (a) v = 0.5 m/s, v = 1 m/s and v = 10 m/s; (b) 
close up of = 0.5 m/s and v = 1 m/s 



give the ship a density of 579 kg/m^. This would give Ps/Pw < 1- However, if the ship was 
filled with water then ps = 1078 kg/m^ and Ps/Pw ^ 1- Also air pockets within the ship 
would effect the sinking kinematics. 

In conclusion, if the Costa Concordia were to suddenly and quickly slip down the steep 
slope she lies on then large waves could be generated and pose a danger to the immediate, 
far field and offshore coastlines. However, given that this is a ship and not a solid mass it is 
likely that the sinking kinematics would not be straight forward. It is difficult to determine 
the density of the ship and how air pockets may effect the way in which she would sink. 
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3. Tsunami inundation 



When it comes to tsunami inundation, one could simply say: rely on the wisdom of 
ancient people. In Japan, for example, hundreds of stone tablets warn citizen about the 
dangers of a tsunami. For example, one tablet in Aneyoshi, a small coastal town, reads: 
"High dwellings are the peace and harmony of our descendants. Remember the calamity 
of the great tsunami. Do not build any homes below this point." As the ancient warnings 
in stone attest, tsunamis are not new to this vulnerable part of Japan. In 1896, during 
the Meiji Period, a tsunami killed at least 22,000 people on the same Sanriku Coast - a 
death toll chillingly close to the recent 2011 disaster which left an estimated 23,000 dead 
or missing. 

3.1. Run-up amplification. Until recently the analysis of long wave run-up on a plane 
beach has been focused on finding its maximum value, failing to capture the existence of 
resonant regimes. Stefanakis et al. [30] performed one-dimensional numerical simula- 
tions in the framework of the NSWE to investigate the Boundary Value Problem (BVP) for 
plane and non-trivial beaches. Monochromatic waves, as well as virtual wave-gage record- 
ings from real tsunami simulations, were used as forcing conditions to the BVP. Resonant 
phenomena between the incident wavelength and the beach slope were found to occur, 
which result in enhanced run-up of non-leading waves (these resonances are different from 
edge- wave resonances and from bay resonances, which are two-dimensional phenomena). 
Run-up amplification occurs for both leading elevation and depression waves. Figure 15 
shows the maximum run-up amplification ratio as a function of non-dimensional frequency 
and non-dimensional wavelength for two beach lengths, namely L = 12.5 m and 4000 m. 
Here rjo is the incoming wave amplitude, uj the frequency, g the acceleration due to grav- 
ity, 9 the beach slope, Aq the incoming wavelength. The classical formula based on linear 
theory for a constant slope reads (see for example [18]) 



Even though these results are in an ideal setting (plane beach, incoming sinusoidal wave) 
they have been shown to persist in more realistic situations (arbitrary bathymetry, complex 
wave signal). The eflFects of dispersion have been studied as well: dispersion reduces 
slightly the amplification ratio. An experimental confirmation of the theoretical results of 
Stefanakis et al [30] was recently published by Ezersky et al [8]. 

3.2. Whirlpool like effect in tsunamis. We conclude this paper by investigating an 
intriguing phenomenon observed during the 11th March 2011 tsunami in Japan. The 
behaviour of the tsunami in Oarai port was picked up by reporters as it seemed to create 
a whirlpool like eflFect within the bounds of the harbour walls (see Figure 16). 

In an attempt to understand this phenomenon better, the VOLNA code was run on a 
simplified geometry of the harbour, inputting a tsunami from the right hand side of the 
domain. The harbour was given a depth of 5 m at the shore with a gently sloping seafioor 
with slope 1/300. These values are based on rough values from Google Earth, as the depth 
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Runup Amplification vs Angular Frequency 
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Figure 15. Maximum run-up amplification ratio as a function of non- 
dimensional angular frequency (top) and non-dimensional wavelength (bot- 
tom) for two beach lengths, namely L = 12.5 m and 4000 m (from [30]). 

is 5 m just outside the harbour and it drops off to 15 m over approximately 3 km. See 
Figures 17 and 18 for a comparison of the real and simplified geometry. The tsunami was 
modelled by a sine wave with amplitude 1 m and period 10.5 minutes 

T] = -sin(O.Olt), at X = 7000 m. 

An unstructured triangular mesh was generated over an approximate area of 7 x 5.5 
km^. The mesh contained approximately 320 x 10^ elements of order 100 m outside the 
harbour and refined to 5 m within the harbour area as shown in Figure 19. 

The free surface results between 7 and 18 minutes are shown in Figure 20. Once the 
wave hits the first harbour wall, it is refiected and refracted. When the wave hits the 
subsequent boundaries a swirling type motion starts to appear (see Figure 21 for a close 
up of this effect between 13 and 16 minutes). This motion is dramatic for a few minutes, 
however even sometime later a depression in the harbour is maintained (see Figure 22). It 
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Figure 16. Image of Oarai Harbour on a Japanese news channel when 
the tsunami hit. Source: http://www.youtube.coin/watch?v=Qgpljq60g4k 

(screenshot 10/08/11) 




(a) (b) 
Figure 17. Satelhte images of Oarai harbour from Google Maps (a) far (b) close. 

is important to point out that the observed vortices are macro-scale ones since the motion 
is assumed to be non-dissipative. More details can be found in [21]. 

4. Conclusion 

This paper shows that there are still a number of poorly explored topics related to 
tsunami science, especially when it comes to tsunami generation or to tsunami run-up. 
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Figure 18. Simplified bathymetry used for Oarai harbour simulation, the 
tsunami is generated at the right boundary and travels from right to left, all 
other boundaries are wall boundaries. 




Figure 19. Unstructured triangular mesh used for Oarai harbour simula- 
tion. There are approximately 320 x 10^ elements of order 100 m, refined to 
5 m within the harbour area. 
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(d) t 15 min (e) t 16 min (f) t 18 min 



Figure 20. Snapshots of the free surface between 7 and 18 minutes. 
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